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Abstract 

Claims have been made that /o(1370) does not exist. The five primary sets of data requiring its existence 
are refitted with suitable Breit-Wigner amplitudes. Major dispersive effects due to the opening of the 4ir 
threshold are included for the first time; the a — > 4-7T amplitude plays a strong role. Crystal Barrel data 
on pp — > 3n° at rest require /o(1370) signals of at least 32 and 33 standard deviations (<r) in 1 So and 3 Pi 
annihilation respectively. Furthermore, they agree within 5 MeV for mass and width. Data on pp — > 777771"° 
agree and require at least a 19a contribution. This alone is sufficient to demonstrate the existence of /o(f370). 
BES II data for J/$ —+ </>7r + 7r~ contain a visible /o(I370) signal > 8a. In all cases, a resonant phase 
variation is required. Cern-Munich data for rnr elastic scattering are fitted well with the inclusion of some 
mixing between a, /o(1370) and / (1500). Values of for / 2 (1565), p 3 (1690), p 3 (1990) and / 4 (2040) are 
determined. 

PACS: 13.25.Gv, 14.40. Gx, 13.40. Hq 

1 Introduction 

The /o(1370) plays a vital role in the spectroscopy of light J p = + mesons. The role of experiment is to 
identify the resonances required by data and determine their parameters: mass, width and branching ratios. 
Many schemes exist for grouping the observed states into nonets. The + glucball is predicted at ~ 1640 MeV 
[1], but with a sizable error from uncertainty in how to normalise the mass scale. The glueball would make an 
additional SU(3) singlet state, mixing with neighbouring nonet states. It is therefore necessary to scrutinise the 
evidence for each state. 

Klempt [2] and Klempt and Zaitsev [3] question whether there is an identifiable + glueball and hence argue 
against the existence of /o(1370). Ochs has argued against it at conferences for two reasons. He suggests that it 
can be confused with a broad 'background' which might obscure the analysis of the 1000-1500 MeV mass range 
[4]; secondly he argues that it has not been identified definitively in Cern-Munich data [5]. The weakness of 
these arguments is that no attempt has been made to refit the data where /o(1370) has been identified. The 
main objective of the present work is to do just that. 

The earliest evidence for the /o(1370) came from experiments on 7T7r — > KK at the Argonne and Brookhaven 
laboratories in the late 1970's. A peak was observed in the S-wave at 1300 MeV in three sets of data: Cohen 
et al. [6], Pawlicki et al. [7] and Etkin et al [8]. In those days it was called the e(1300). Further evidence 
appeared in the years 1992-6 from several experiments in quick succession. Amsler et al. reported a peak in the 
7777 S-wave at ~ 1400 MeV in Crystal Barrel data on pp annihilation at rest into 777771"° [9] (as well as a further 
peak at 1560 MeV, later identified as /o(1500)). Gaspero reported a + state in 47r at 1386 ± 30 MeV in pd 
bubble chamber data at rest in the reaction pn — > 27r + 37r _ [10]; the Obelix collaboration quickly confirmed this 
observation in rip — > 37r + 27r _ [11], quoting a mass of 1345 ± 12 MeV. Crystal Barrel reported a similar + state 
in 4tt at 1374 ± 38 MeV in pp annihilation at rest to (7r + 7r~7r°7r°)7r [12]. 

The earliest fits in which /o(1370) and /o(1500) appeared together were made using early low statistics 
Crystal Barrel data on pp — ► 37r° at rest [13], [14], [15]. Publications using the full statistics data (used here) 
appeared in 1995 [16] and 1996 [17]. The latter work involved a simultaneous fit [18] to Cern-Munich data on 
7T7T — * 7T7T [19] and data on 7T7t — > KK. The observation of nearby /o(1370) and /o(1500) states, both dominantly 
non-strange, excited interest in the + glueball. Further extensive studies of /o(1370) and /o(1500) by Crystal 
Barrel, Obelix and Omega collaborations may be traced via the Particle Data Tables [20]. 

The present work fits simultaneously the 5 most definitive sets of data available to me. The first two are 
Crystal Barrel data on pp — > 37r° at rest in liquid hydrogen [17] and gaseous hydrogen. These two sets of data 
allow a clean separation of annihilation from 1 So and 3 Pi pp states: P-state annihilation is ~ 13% in liquid and 
48% in gas. The combined fit also includes two sets of data for pp — > 777771" in liquid and gas. There are visible 
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/o(1370) and / (1500) peaks in both of these sets of data. The fifth definitive set of data comes from BES II for 
J/\]/ — > (f>n + ir~ . There is a visible tttt peak at 1350 MeV, attributed to interference between / (1370), /o(1500) 
and /2(1270) [21]. Those data are refitted with and without inclusion of /o(1370). 

Finally, an earlier analysis fitted data on tttt — > KK and r\r\ [22]; this analysis also fitted Kloe data on 
(f) — > 7r°7r°7 and t]tt°j. The tttt — ► KK data suggest the presence of /o(1370) but cannot be considered definitive 
compared with the other sets of data. This analysis remains consistent with parameters of /o(1370) found here. 

The Particle Data Group quotes very large errors for mass and width of /o(1370): (1200— 1500) — i(150— 250) 
MeV. From the present analysis, the error on the mass is very small: ±15 MeV from systematics. The resonance 
will still be referred to as /o(1370), despite the fact that its peak position in the tttt channel comes out nearly 100 
MeV lower. The Particle Data Group appears to be influenced by large variations in masses and widths fitted to 
Att data. These variations are due to the fact that there has been no serious attempt to include a — > Att in these 
analyses. A major objective here is to treat fully the dispersive effects due to the opening of the Att threshold, 
to which a, /o(1370) and /o(1500) all couple strongly. 

The full form of the Brcit-Wigner resonance formula, 



We shall find that m(s) is large, indeed larger than (M 2 — s) in the Breit-Wigner denominator. It turns out that 
this places severe limitations on the ratio T2t T /T4 7T which can be fitted to data. This point was not appreciated 
in earlier work. Nevertheless, a good solution emerges naturally. 

Here there is an interesting conclusion. The fit includes explicitly the s-dependence of the Att channel, as far 
as present data allow. However, the amplitude still produces an Argand diagram very close to a circle and hence 
similar to a simple pole. This justifies to some degree the common usage of a simple Breit-Wigner formula in 
fitting data. 

A third point concerns the coupling to Att by the broad component in the tttt S-wave related to the a pole; 
this component will be called a as a short-hand. It plays an essential role in fitting all data where a 7T7t pair is 
produced. Let us review the situation briefly. 

The tttt S-wave phase shift reaches 90° at <~ 900 MeV. There is an Adler zero in the clastic amplitude 
at s ~ m2 -/2, just below the tttt threshold. The resulting amplitude rises approximately linearly with s. In 
production data from E791 [24] and BES II [25], a strong peak is observed at <~ 500 MeV. Both these production 
data and elastic scattering may be fitted with the same Breit-Wigner denominator, but with an s-dependent 
width and with different numerators for production and elastic cases; for elastic scattering, the numerator 
contains the Adler zero, whereas for the BES II production data fitted here it is consistent with a constant. The 
same variation of the phase shift with s is observed in elastic scattering and production data from 450 to 950 
MeV within errors of ~ 3.5° [26]. The large displacement of the pole from 900 MeV arises from the s-dependence 
of the width; the Cauchy-Riemann relations control the s-dependence of the real and imaginary parts of the 
amplitude as one extrapolates from the physical s-axis. 

The fit of Ref. [22] to tttt — ► KK and r]rj data and the Kloe branching ratio for <p — * 7r°7r°7 determines ratios 
of coupling constants for the a coupling to KK and 7777: g 2 (KK)/g^ = 0.6±0.1 and g^/g% v — 0.2±0.05. The 
remaining unknown is the coupling of a to Att. This will play an essential role in the work reported here. 

A vital question is whether this coupling to Att eliminates the requirement for the /o(1370). Tornqvist has 
suggested [27] [28] that a second pole could appear in the a amplitude due to the opening of the Att threshold. 
Could this explain /o(1370) as a non-gg state? A related point is that Maiani et al. have suggested that /o(1370) 
may be a molecular state [29]. 

The layout of the paper is as follows. Section 2 discusses the formalism and gives equations. Amplitudes are 
expressed in terms of T-matrices for reasons discussed there. This raises some issues concerning how to fit elastic 
scattering. Readers interested only in results may skip this section, but the issues going into the formulae are 
outlined in words for the general reader. It is important to add that extensive fits to Crystal Barrel and other 
data have been made by Anisovich and Sarantsev, using K-matrix techniques [30]. These analyses produce results 
closely similar to the present work, and there is no conflict between the two analyses: they should be regarded as 
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contains a real dispersive term m(s) [23], which for the Att channel reads 
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complementary views. The K-matrix analyses use several sets of data not available to me, for example Crystal 
Barrel data on pp — > KKir. Their conclusion is that /o(1370) is needed, with mass and width in close agreement 
to what is found here. However, they do not address the question of how much the fit changes if / (1370) is 
omitted. 

Section 3 is the heart of the paper, concerning fits to data with and without /o(1370). A suggestion made by 
Ochs is that the amplitude for the tttt S-wave should be fitted freely in magnitude and phase in bins of tttt mass, 
without assuming a Breit-Wigner form. Over the limited mass range 1100-1460 MeV, this is done and confirms 
the assumption of a resonance form for the amplitude. 

Section 4 describes the simultaneous fit made to Cern-Munich data on tttt elastic scattering. They can be 
fitted slightly better with / (1370) than without, but cannot be considered definitive on this question. Section 
5 describes the fit to 777771-° data; there is a visible peak in these data, sufficient alone to justify the existence of 
/o(1370). 

Section 6 describes the fit to BES data on J/* — > </>7r + 7r~. An important detail is that these data, together 
with data on J/* — > uitttt [25] and loKK [31] require the existence of an / (1790) distinct from / (1710). The 4>tttt 
data contain a clear tttt peak at 1790 MeV. Data on <f>irir and tpKK final states require a ratio £?i?[/o(1790) — > 
tttt}/ BR[j "0(1790) -> KK] > 3 [21]. Data on J/* -> loKK contain a strong KK peak due to / (1710); this peak 
is completely absent from ujtttt data which require Bi?[/ (1710) -> 7r7r]/B J R[/ (1710) -> KK] < 0.11 with 95% 
confidence. A single resonance must have the same branching ratios in all sets of data, whereas these data differ 
in branching ratio by at least a factor 22. So two separate resonances are required. Unfortunately, the PDG 
ignores these decisive results and continues to lump / (1790) with / (1710). In the present fits to pp — > 3ir° 
data, /o(1790) is included together with / (1370) and /o(1500), though its eventual contribution to the data is 
small. 

Section 7 makes brief remarks on data for tttt — > KK, and section 8 remarks on the need for further analyses 
where the opening of thresholds leads to large dispersive effects. Section 9 makes concluding remarks. 

2 Discussion of Formulae 

There are two obvious ways of expressing amplitudes, using iiT -matrices or T-matrices. Each has some advantages, 
but also limitations. It seems likely that neither is perfect, so approximations are needed in cither case. 

Resonances appear directly as poles of the T- matrix. The tttt S-wave is of primary concern. For this amplitude, 
if-matrix poles are displaced strongly from T-matrix poles. For elastic scattering, K cx tan 8 (where 8 is the 
phase shift), and if-matrix poles are at <~ 700 and 1200 MeV, whereas the /o(980) pole is at 998 — H7 MeV. 

In fitting production data, e.g. pp — > 3tt°, it is obviously advantageous to use T-matrices, because resonances 
are T-matrix poles. The primary objective of the present work is to test whether the /o(1370) is needed or 
not. It is necessary to move its mass, width and couplings to all decay channels in steps, so as to examine the 
effect on x 2 - This cannot be done readily using if-matrices, since the /o(1370) is built out of a combination of 
if-matrix poles. It is also necessary to remove /o(1370) from the fit and again test \ 2 . This cannot be done in 
a controlled way using ^-matrix poles: if one if-matrix pole is removed, all resonances are affected. For this 
reason, amplitudes will be written directly in terms of T-matrices. 

There is a second related point. It is well known that minimisation routines converge best when expressed 
in terms of eigenvectors, i.e. T-matrix poles. Weak or questionable resonances appear as weak eigenvectors and 
can be recognised immediately from the error matrix of the fit. 

Yet another consideration is that earlier work fitting the a pole and Cern-Munich data was done using T- 
matrices [32]. The formulae used there are readily expanded to incorporate the Att channel. Also Ref. [22] fitted 
data on tttt — > KK and 7777 and also Kloe data with T-matrices; it is valuable to maintain consistency with that 
analysis, for comparison of results. 

2.1 Elastic Scattering 

There are however questions about how to deal with elastic scattering. Below the KK threshold, the amplitude 
is confined to the unitary circle. Both /o(980) and a contribute, as do the low mass tails of further resonances. 
Experiment shows directly how to treat the overlap of these resonances. In Cern-Munich data, the tttt phase 
shift rises dramatically near 1 GeV from <~ 90 to 270° due to the narrow /o(980). The appropriate treatment 
below the KK threshold is to add phases, hence multiply S'-matrices: S — exp(2i<5). 
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Above the inelastic threshold, multiplying S'-matrices gives a fit of rather indifferent quality. It is clear 
that other factors must be relevant. If one solves a relativistic Schrodinger equation using a trial potential which 
reproduces a, /o(980) and /o(1370), the solution is automatically unitary and analytic. One finds that amplitudes 
differ from both (i) the product of individual T-matrices for each resonance, (ii) the sum of if -matrices. Neither 
gives an accurate parametrisation. The main problem appears to be that resonances mix via processes of the 
form < cr|7T7r|/o > or other intermediate channels KK , 7777, etc. Mixing is strong for clastic scattering, since the 
amplitudes are at the unitary limit when one takes account of all channels. The mixing gives rise to well known 
level-repulsion. This repulsion is still highly significant one full-width away from the resonance mass. 

Formulae for mixing have been given by Anisovich, Anisovich and Sarantsev [33] and will be reproduced here 
in a slightly modified notation. For the 2-resonance case, the denominator D(s) of one resonance may be written 

DuM = Ml - . - mi (») - .M.rrM - w -S^M- m . (3) 

Mixing arises via £12 which in general is complex and may be s-dependent. The propagator matrix describing 
two resonances is then 
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Because a overlaps strongly with /o(980), /o(1370) and /o(1500), this mixing has been included explicitly in 
fitting elastic data. It turns out that inclusion of mixing between these pairs of states leads to an excellent fit 
using constant values of Bi 2 . 

It is also instructive to expand the denominator of the last term of eqn. (3) off resonance using the binomial 
theorem. The result is 

£ 12 B 21 /. m 2 (s)+iM 2 r*°* > 
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Mf 



Mf 
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From the last term, one gets contributions of the same sign to ReDn(s) and ImDn(s), while from the first term 
[Mf — s — imi — iMiT\ ot ] contributions have opposite signs. The result is to rotate the phase of the amplitude, 
which derives purely from D(s). This rotation is large when resonances overlap strongly. Unless the mixing is 
included explicitly, one must expect that the resonance denominator may need to be multiplied by a phase factor 
exp(i^). This was indeed observed in Ref. [22], where it was sufficient to take <p as constant. This point will be 
relevant in fitting 1~~ and 2 ++ states to elastic data. 



2.2 The continuation of amplitudes below thresholds 

Consider as an example tttt — > KK. It is necessary to continue this amplitude below threshold; in T-matrix 
language this is the analytic continuation of the T i2 component. As a result, there are contributions to tttt elastic 
scattering from sub-threshold tttt — > KK — > tttt. However, caution is needed in making this continuation. The 
tttt — ► KK amplitude is proportional to phase space p 2 = yl - 4M|-/s for the KK channel. Below threshold, 
this continues analytically as iy/AM^/s — 1. The analytically continued amplitude rises rather strongly below 
threshold. If one is not careful, this continued amplitude can make dominant contributions to elastic scattering 
below the inelastic threshold. This is counter-intuitive. A particular case arises for / 2 (1565), which couples 
strongly to luuj. If one does nothing about the large sub-threshold ww contribution, it produces big interferences 
wth the nearby / 2 (1270) and can distort rather severely the mass and width fitted to / 2 (1270). 

The answer to this point is straightforward. Amplitudes above threshold contain form factors due to the 
finite radius of interaction forming the resonance. In Ref. [22], data on tttt — > KK were fitted empirically to an 
exponential form factor exp(— 5.2k 2 ), where k is momentum in the KK rest frame in GeV/c; a study of / 2 (1565) 
in pp — > (louj)tt at rest in Crystal Barrel data also requires a form factor above the uiu> threshold of similar 
strength [34]. The analytic continuation can be evaluated below threshold using a dispersion integral. This is 
not an accurate procedure because of uncertainties in the tttt — > KK amplitude above the available mass range. 
However, the qualitative feature emerges that a rapidly falling form factor is required below threshold as well as 
above. In Ref. [22] this sub-threshold form factor was fitted to Kloe data with the result cxp(— 8.4|fc| 2 ). This 
empirical cut-off will be adopted here for sub-threshold coupling to KK, 7777 and uu. 
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Figure 1: At: phase space for pp (dashed), era (chain curve) and the fit adopted here (full curve) 



2.3 Treatment of the 4-7T channel 



The An phase space volume may be modelled [18] by the production of two resonances ( pp or acr): 

I l,n- 1" Jim 2 K y/SSlS 2 



P4tt(s) = / 

J4r, 



inis^T^fFFis), 



(6) 



where p\ and p 2 are momenta of pions from decays of each resonance in its rest frame, and p stands for the 
momenta of the p or a in the centre of mass frame. In Ref. [23], extensive illustrations are shown of the 
dependence of p4?r(s) on exponential form factors FF. These factors begin to play a significant role at ~ 1.45 
GeV. A form factor 

FF = cxp[-(s- 1.45 2 )] (7) 

is chosen in present work with s in GeV 2 . If one were fitting data on nn — > Air, this form factor would be 
rather important. However, for the present study of /o(1370) and /o(1500) — ► tctc it has only rather small effects 
within errors. This is because the term in the Breit-Wigner denominator cuts off the tttt channel strongly 
at high mass; the /o(1370) — ► -kit amplitude is already quite small at 1.45 GeV, where the form factor begins. 
The /o(1500) is sufficiently narrow that the effect on the line-shape from An inelasticity is rather small. The one 
place where the form factor is important is in suppressing high mass contributions to m(s). 

Fig. 1 shows pp phase space as the dashed curve and era phase space as the chain curve. They are rather 
similar. Their relative contributions to each resonance are poorly known. The strategy here is to parametrise 
P4ir empirically as 

Pi "~ l + exp[A( So -s)r ( ) 

The parameters A and so in the Fermi function of the denominator are optimised in the overall fit, with the 
result A = 3.39 GeV -2 , s = 3.238 GeV 2 . The result is shown by the full curve on Fig. 1. 

The dispersive contribution m(s) to Breit-Wigner amplitudes is evaluated numerically at 10 MeV steps of 
mass, and the programme interpolates quadratically in mass using the nearest 3 bins. The dispersion integral 
is the same for a and all /o's, except for (i) a subtraction at the resonance mass M, where the real part of the 
amplitude is zero, (ii) a scaling factor depending on the coupling constant to 4w. The sub-routine for doing the 
principal-value integral is available from the author if it is needed in other cases. 



2.4 Explicit equations for resonances 

Formulae for the a amplitude follow the same form as used in earlier work on the a pole [32], except for the 
inclusion of m(s) for Att. Equations will be repeated here for completeness. The elastic amplitude is written 

Tn (a) =N{s)/D{s). (9) 

The numerator contains an Adler zero at s — sa — 0.41m 2 . For fits to BES data on J /ty — > uirir, the numerator 
is taken as a constant. Channels 7T7t, KK, -q-q and At: will be labelled 1 to 4. The propagator of the a is given by 

D(s) = M 2 - s - g\ S ~ SA z s - m(s) - iMT tot (s) (10) 
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MTi{s) = gl i^- A pi{s) 

g \ = M(b 1 +b 2 s)exp[-{s- M 2 )/A] (12) 



3i(s) = 



2 + p\ln e 



I -Pi 



(13) 



1 + 

z s = h{s)-h{M 2 ) (14) 

MT 2 {s) = 0.6g 2 FF 2 (s) (15) 

MT 3 (s) = 0.19g 2 FF 2 (s) (16) 

F^( S ) = exp(-a|fc|, 2 ) (17) 

Mr 4 (s) = Mg iP ^{s)/p^(M 2 ). (18) 

The value of a is 5.2 above thresholds and 8.4 below. The numerical coefficient 0.19 in eqn. (16) has been revised 
very slightly using the branching ratio fitted in Section 5 between rjrj and irir. Values of b\, b 2 , A and M of eqn. 
(12) are given below in table 6. 

Resonance denominators for /o(1370), /o(1500) and /o(1790) are taken in the form of Eq. (1). An important 
detail is that the factor (s - s A )/(M 2 - s A ) of Eq. (10) is also used for Ti, T 2 and T 3 of / (980), / (1370), 
/o(1500) and /o(1790). That is, the Adler zero is included into the widths of all + resonances. Otherwise, 
parameters of /o(980) are taken initially from the BES determination [21], but are re-optimised within the sum 
of statistical and systematic errors when fitting present data. For /o(1500), the ratio r^/Ti is taken from the 
PDG average. For /o(1370), it is taken from Ref. [30], where extensive Crystal Barrel data on pp — > KKtt are 
fitted. However, in practice the ratio T 2 /Ti has rather have little effect here; only the full width of the resonance 
is crucial. Values of r 3 /ri are determined in Section 5, but again have little effect. The same ratios are used for 
/o(1790) as for / (1500), in the absence of good data for / (1790) in the KK channel; since the /o(1790) signal 
is weak, this is of no consequence. A trial has been made including into /o(1500) a weak coupling to lolo with a 
coupling constant a third of that for known decays to pp; the effect is negligible. 

The /2(1270) is parametrised using 

r -<»» - r -< M2 »Sai < 19 > 

r„M = r„(M 2 )-«^i-. (20) 

Here k is the pion momentum in the f 2 rest frame and k r is the value on resonance. The D 2 are Blatt-Weisskopf 
centrifugal barrier factors. Expressions for them are given in Ref. [18] at the end of Section 2.1; the barrier 
radius optimises at 0.75 ± 0.04 fm. The KK and r/r) channels are treated in the same way. From the Particle 
Data book, r 27r /r to t = 0.847 on resonance. This value is not sufficiently accurate for fitting Cern-Munich data 
and is re-optimised to 0.802, since the relative heights of p(770) and /2(1270) are important. This value may be 
accounting for a form factor in irp — > mrp over the mass range between these two states. 

The p(770) is parametrised like /2(1270) using the D\ centrifugal barrier factor. Its mass and width optimise 
at 778 and 153 MeV, close to PDG values. Trials were made including coupling to KK and W7r, but had no 
significant effect on the fit to Cern-Munich data compared with uncertainties in the p(1450) contribution. 

The /2(1565) presents a problem. It is known from Crystal Barrel data [34] that it couples strongly to lolo. 
It is expected to couple to pp a factor ~ 3 more strongly. However, there are unfortunately no data to confirm 
this strong pp coupling. Consequently it is dangerous to fit /2(1565) with the full s-dcpcndcnce and m{s) used 
for /o's. In view of the fact the /o(1370) and /o(1500) line-shapes come out close to those of simple poles, it is 
therefore fitted with a Brcit-Wigner amplitude of constant width, except that the lolo channel is added explicitly 
This channel has an important effect in cutting off the line-shape above resonance, see Fig. 9 below. 

Other 1~~, 2 ++ , 3~~ and 4 ++ resonances are likewise fitted with Breit-Wigner amplitudes of constant width, 
for simplicity and ease of comparison with other work. In fitting elastic data, each amplitude is multiplied by 
a fitted factor exp(icj>). This is done rather than including mixing between all resonances, since the number of 
mixing terms becomes too large, and furthermore the resonances are poorly separated at present. 
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3 Fit to pp — > Stt at rest 

Fig. 2 shows the Dalitz plot for data in liquid hydrogen at rest; it is similar for gaseous hydrogen. [This figure 
is from Ref. [16] where /2(1520) = /2(1565)]. Mass projections for both sets of data are shown in Fig. 3. There 
are conspicuous peaks due to /2(1270) and /o(150Q) + /2(1565). Narrow bands due to / (980) are also visible. 
Ingredients in the fit are a, / (980), / 2 (1270), / 2 (1565), /o(1370), / (1500) and a weak / (1790). 

P-state annihilation is best determined by data in gas, where it makes up 48% of events, close to the predicted 
50% from calculations of Stark mixing [35]. This is mostly annihilation with orbital angular momentum L = 1 
in the transition, but a minor detail is that there is a little L = 3 production of 7r°/2(1270) and 7r°/2(1565). The 
fit to liquid data uses 3 Pi and 3 P2 components scaled down between gas and liquid by a factor which optimises 
at 0.155. The eventual P-state fraction is 12.8% in liquid. 

A warning is that amplitudes for 1 So, 3 Pi and 3 P2 — > /2(1270) or / 2 (1565) are poorly separated without 
using information from interferences between the three 7r°7r° contributions, which will be labelled 12, 23 and 13. 
The full angular dependence of amplitudes involves d-matrices which take account of rotations of axes between 
12, 23 and 13 combinations: 

3 Fi,J z = +l : 3 cos a sin 6> cos 6> - sina(3cos 2 9 - 1) (21) 

3 Pi,J 2 = : V2[3 sin a sin 9 cos 9 + cosa(3cos 2 0—1)] (22) 

3 P2,Jz — 2 '■ cos a sin 2 9 — sin a sin 9 cos 9 (23) 

3 P2,Jz — l '■ sinasin 2 9 + cos a sin cos 0. (24) 

For the 12 combination, 9 is the angle in the 7r°7r° rest frame between 7Ti and the recoil 7^; ai_3 are lab angles 
of 12, 23 and 13 combinations in their plane. 

It turns out that 3 P 2 annihilation dominates strongly over 3 Pi for both /2(1270) and / 2 (1565). Visible 
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evidence for 3 P 2 - * /2(1565) is the enhancement in Fig. 2 near the centre of each band just above 1500 MeV; 
this is how the Asterix collaboration discovered /2(1565) in gas data [36]. If 3 Pi annihilation to / 2 (1270) were 
large, there would be strong constructive interference between any two crossing bands near cos 9 = 0.6; there is 
no sign of any such enhancement in the data. 

3.1 The a amplitude in 3ir° data 

The a is very broad, extending over the full tttt mass range from 0.27 to 1.74 GeV. Over such a large range, some 
s-dependence is to be expected in the numerator N(s) of the production amplitude. The s-dependence required 
by the data is quite different between elastic scattering and the production process. An important feature of 
the data is an area of low intensity at the centre of the Dalitz plot of Fig. 2. Fitting this feature is delicate and 
demands s-dependence in N(s). 

Extensive trials have been made using a numerator for the a amplitude of the form 

N(s)=A[l + Bs + C/(s + s Q )], (25) 

with so > 0. What emerges is that (a) either complex B or complex C is definitely required and can fit the 
data well, (b) including both B and C over-parametrises the amplitude, i.e. strong correlations develop between 
B, C and s . It is better to tolerate a small increase in x 2 , so as to keep essential features clear with minimal 
correlations between fitting parameters. Final fits were made with the form A{1 + Bs). Tests were also made 
including amplitudes arising from the opening of the KK threshold, i.e. oc T nK . These turned out to be 
negligible. A fit was also tried using a dependence on spectator momentum k given by the diffraction pattern of 
a black disk, but this gave a poor fit. 

3.2 Prelimiminary remarks on the goodness of fit 

The data are available to me only in the form of binned data, rather than individual events. Initial fits revealed 
that some edge bins have abnormally high % 2 . All lie immediately at the edge of the Dalitz plot. This is actually 
visible on Fig. 2. The same effect is observed for 777777° data. The obvious explanation of these bad bins is that 
the acceptance near the edge of the Dalitz plot may be incorrectly assessed. These bad bins have been removed 
without any significant effect on fitted amplitudes. 

The resulting \ 2 is 2.83 per bin for 3tt° data in liquid, 2.85 in gas, 2.64 for 777777° in liquid and 2.70 in gas. A 
close inspection of discrepancies over the Dalitz plots reveals only an apparently random scatter, with no clear 
systematic effects. The mean x 2 is 2.85 per bin. To allow for this, all values of % 2 quoted in the paper are scaled 
down by a factor 2.85 so that the average % 2 becomes 1 per bin. This is necessary for a correct assessment of 
the significance level of observations. 

3.3 Fits including / (1370) 

The Argand diagram for one 7r°7r° 1 So combination is shown in Fig. 4(a); individual a and /o(1370) components 
and their coherent sum are shown in other panels. At low mass, there is a conspicuous loop which is well fitted 
with the a pole. This feature was correctly diagnosed by Ishida et al. [37]. They fitted only the 7T7T mass 
projection, so the present fit is much more accurate. 

At higher mass, there is a loop due to / (980). Here a detail needs explanation. Fig. 5 shows the Argand 
diagram for / (980) alone in tttt — ► tttt. A mean kaon mass of 495.7 MeV is assumed. The KK inelasticity sets 
in very rapidly at threshold, and the peak inelasticity is at 1.010 GeV. Thereafter, the inelasticity parameter 77 
rises again slowly, The result is a definite 'dent' in the Argand diagam at 1.01 GeV. 

Returning to Fig. 4(a), there is a further loop at <~ 1300 MeV, followed by a large loop due to / (1500). 
Fig. 4(a) resembles closely fits made in 1996 [17]. The loop at 1300 MeV is the feature which is crucial to the 
existence of /o(1370). The vital questions are: 

• (i) is this loop really needed? 

• (b) could it be fitted with a and /o(1500) alone, without the need for / (1370)? 
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Figure 4: Argand diagrams for the irir S-wave in liquid hydrogen; masses are marked in GcV 
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Figure 5: The Argand plot for / (980) alone; masses are marked in GeV 



3.4 A suggestion of Ochs 

Ochs has questioned whether the loop at 1300 MeV is an artificial feature of the way amplitudes are parametrised 
in terms of resonances (or backgrounds). My comment is that the parametrisation must be analytic; resonance 
forms used here satisfy this condition. Nonetheless, he has suggested fitting the irir S-wave amplitude freely in 
magnitude and phase in bins of tttt mass, to see how definitively the data require the loop at 1300 MeV. 

It is not possible to do this over the entire Dalitz plot, because of strong interferences between one low mass 
7T7T contribution and two at higher mass. It is necessary to rely on the strong a loop at low mass and also the 
existence and parameters of the /o(1500), which is well known today from other data; it is also important to 
constrain the /o(980) within bounds set by BES II data. 

It is, however, straightforward to fit the S-wave amplitude freely in magnitude and phase in bins from 1100 
to 1460 MeV, i.e. over the 1300 MeV mass range. In this test, parameters of /o(1500) and /o(980) are allowed to 
re-optimise within the narrow ranges allowed by other data. For the a, coefficients A and b of N(s) are set free, 
but hardly move because they are determined by elastic data. Results are shown in Fig. 6 by crosses (indicating 
errors). Deviations from the fit of Fig. 4(a), shown by the full curve, are barely above statistics, except for the 
third point (labelled 3) which moves by ~ 2.5 standard deviations. The agreement between the binned fit and 
Fig. 4(a) rules out the possibility that the 1300 MeV loop is an artefact of the parametrisation. As one example, 
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Figure 6: The Argand diagram for the mr S-wave (full curve) compared with free fits to magnitude and phase 
in nine 40 MeV wide bins of mass from 1100 to 1460 MeV (crosses) 



it is not possible to replace the 1300 MeV loop by a narrow cusp arising from interference effects; such a cusp 
will be illustrated below for the mr D-wave. 



3.5 Properties of the fitted / (1370) 

The fitted dispersive contribution m(s) for /o(1370) is shown in Fig. 7. Near 1300 MeV, it varies roughly linearly 
with s, like the term (M 2 — s) in the Breit-Wigner denominator. However, m(s) is a factor ~ 1.6 larger than 
(M 2 — s), This in an unusual feature, showing that care is needed in treating the 4ir threshold correctly. It also 
leads to some effects which had not been forseen. 




1 .5 2 
Mass(GeV) 

Figure 7: The dispersive term m(s) for / (1370) 



The term m(s) is directly linked to MT^(s) by the dispersion relation eqn. (2) of Section 1. This relation 
is reproducing effects of the loop diagram of Fig. 8. Results are analogous to vacuum polarisation, and lead to 
renormalisation effects in the Breit-Wigner denominator if the magnitude of is scaled up or down. 



Figure 8: The dispersive term m(s) for / (1370) 



A consequence in the combined fit to pp — ^ 3tt° and elastic data is that there is a tight constraint on the ratio 
^2tt/^4:tt- Conversely, there is considerable flexibility in the absolute value of r 27r which can be fitted. Table 1 
shows pairs of values fitted to data; T 27r and are almost linearly related until r 27r approaches zero, when the 
fit deteriorates rapidly. The final fit uses the lowest value of T 2 n giving a satisfactory fit, namely 325 MeV. A 
surprise is that /o(1370) is fairly elastic on resonance, though the inelasticity increases rapidly thereafter. Over 
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the range of values shown in Table 1, there is almost no visible change in the line-shape of / (1370). This is a 
renormalisation effect. It is illustrated in Fig. 9(b) for two widely different values of r 27r . 
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Tabic 1: Parameter variations with T 27I of / (1370); units are GeV. 
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Figure 9: Line-shapes of a, /o(1370), /o(1500) and /2(1565) in nn — > tttt, normalised to 1 at their peaks; the 
dashed curve for /o(1370) shows the effect of changing T2-K from 325 to 800 MeV 



3.6 Lineshapes 

Fig. 9 shows the line-shapes of a, /o(1370), /o(1500) and /2(1565) in tttt elastic scattering. The /o(1370) is 
almost degenerate with / 2 (1270), which explains why it has been hard to observe experimentally. The peak of 
the intensity is at 1282 MeV in elastic scattering, though the phase shift goes through 90° at M = 1309 MeV. A 
similar pole mass is quoted by Anisovich and Sarantsev [30]: 1306 — il47 MeV for their solution 1 and marginally 
different for two alternative solutions. Although their parametrisation of amplitudes in terms of if-matrices is 
quite different to the approach adopted here, it is clear that their eventual fit is very similar to mine, and there 
is no serious disagreement between the two types of formalism. 

The /o(1370) and /o(1500) are asymmetric, because the rapidly rising in the Breit-Wigner denominator 
cuts off the intensity at high mass. The /2(1565) is likewise cut off strongly at high mass by the opening of the 
lulu threshold. It peaks exactly at the lulu threshold and has a half-width of 131 MeV below resonance. 

3.7 The 7777 D-wave 

Fig. 10 shows the Argand diagram for 1 Sq annihilation to the 7T7T D-wave. An interesting feature, already noted 
in Ref. [17], is the appearance of a cusp at 1420 MeV, almost midway between / 2 (1270) and j2(1565). No 
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Figure 10: The Argand plot for the irir D-wave in 3ir° data in liquid hydrogen; masses are shown in GeV 



resonance can be accomodated at 1420 MeV by fitted magnitudes and phases of the amplitude. However, the 
PDG lists a state /2(1430). It seems likely that this is an artefact due to similar cusps in other reactions; it is 
an acute observation on the part of the experimental groups finding the effect. 

3.8 P-state annihilation 

Fig. 11 shows Argand diagrams for 3 Pi annihilation to one irir S-wave combination. The result comes almost 
purely from data in gas. There is some similarity to results in liquid, but definite differences; for example, 
the /o(980) contribution is almost negligible. There is again a loop near 1300 MeV in Fig. 11(b). However, 
interference between / (1370) and /o(1500) plays a decisive role in obtaining a good fit. 
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Figure 11: Argand diagrams for the irir S-wave in 3 Pi; annihilation; masses are shown in GeV 



3.9 Fits without / (1370) 

The entire fit has been re-optimised without /o(1370), refitting all parameters. The result is that (renormaliscd) 
X 2 increases for 3ir° data in liquid by 1040, and by 1088 for data in gas. So the significance level of / (1370) is 
slightly over 32 standard deviations in liquid and close to 33 in gas. These are essentially independent, since 1 So 
annihilation to /o(1370) is determined almost purely by liquid data and 3 Pi annihilation by gas data. 

Attempts were made to improve the fit using different parametrisations for the a, e.g. N(s) cx A[l+C/(s+so)]\ 
however, that choice actually gave a slightly worse fit and other types of fit gave no significant improvement. 

Fig. 12 shows the resulting Argand diagrams for the irir S-wave in (a) liquid, (b) gas. The fit tries to remedy 
the situation by using the a amplitude as a replacement for / (1370); the main latitude lies in increasing the 
inelasticity of a to 4ir. This fails because the amplitude cannot move round the required loop quickly enough, 
resulting in a very poor x 2 ; one can see on Fig. 12 that the size of the loop at 1300 MeV has increased compared 
with that of Fig. 4. For the a amplitude, the inelasticity to 4ir rises slowly over the entire mass range 1200 to 
2100 MeV. This produces a slow loop on the Argand diagram, see Fig. 4(b) for liquid(L). It is not possible to 
fit the narrow / (1370), whose full-width at half-maximum is 207 MeV, with this slow loop. 
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Figure 12: Argand plots for the mr S-wave in 3tt° data in liquid (L) and gas (G) with / (1370) omitted; masses 
are shown in GeV 



Tornqvist [?] has remarked that the dispersive contribution m(s) to the a amplitude could induce a second 
a pole near the 47r threshold. All fits have been examined for such a pole, but there is no trace of it in any fit. 

The conclusion is that a narrow /o(1370) is highly significant. However, one should not rely purely on x 2 . 
What adds considerable confidence is that fitted values of mass and width are in excellent agreement between 
two almost independent determinations in liquid data and gas. The value of r 27r is held fixed at 325 MeV in 
both cases. Then the fitted value of on resonance changes by only 4 MeV. Parameters are shown in Table 2. 
One must add a systematic error common to both determinations; the systematic errors shown in Table 2 cover 
the entire range of all observed fits to the six sets of data with any parametrisation of N(s) for the a amplitude. 
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1165 and 1372 


1418 and 1540 
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Table 2: Parameters in MeV for /o(1370) in (a) liquid, (b) gas, (c) combined, and /o(1500) for the combined fit. 

For those who wish to reproduce /o(1370) and /o(1500) with less elaborate formulae than used here, the 
advice is to aim to reproduce the peak and half- heights of Table 2. 

A further test has been made removing the phase variation of /o(1370). There is an increase in x 2 of 165, 
i.e. nearly 13 standard deviations. However, the magnitude of the fitted /o(1370) increases by a large amount 
and the a component decreases. So the non-resonant /o(1370) is then obviously simulating a large part of the a 
amplitude. This confuses the interpretation of this test. 

3.10 / (1370) line-shape in 4tt 

Fig. 13(a) shows as the full curve the line-shape of /o(1370) in mr — » mr; the dotted curve shows the line-shape 
of a simple Breit-Wigner amplitude of constant width with the same mass and width on resonance. The dashed 
curves show estimates of what is predicted in mr — > 4ir. There is some uncertainty concerning form factors at 
high mass. Fortunately these form factors play little role in fitting present data, because the amplitude in the 
2tt channel is already quite small at 1.45 GeV, where the form factor begins to have an effect. The three curves 
illustrate results using form factors exp— a(s — 1.45 2 ) (with s in GeV 2 ); the top curve is for a = 0.5, the central 
one for a = 1.0 and the bottom one for a = 1.5. These span a reasonable range of possibilities. 

The main conclusion from these curves is that the peak of the resonance moves quite strongly between 27T 
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Figure 13: (a) line-shapes of /o(1370) for 2n (full curve), a Breit-Wigner amplitude with constant width (dotted), 
and for An (dashed), (b) the phase angle measured from the bottom of the Argand plot (full curve) and for a 
Breit-Wigner amplitude of constant width (dashed); horizontal lines mark phase shifts of n/2 and n, (c) Argand 
plot; masses are shown in GeV. 



and An channels because of the difference between p^ and p^. The 2n peak is at 1282 MeV and that in An is 
at 1331 MeV, i.e. a mass difference of ~ 50 MeV. However, if one works from the half-heights of the peaks, the 
difference is larger. For the 2n channel, half-height is at 1165 and 1372 MeV, i.e. a mean of 1269 MeV and a 
full-width at half maximum (FWHM) of 207 MeV. For the Air channel, the corresponding half heights are 1241 
and 1514 MeV, i.e. a mean of 1377 MeV and FWHM = 273 MeV. Further data on nn -» An would be very 
valuable. In particular, data separating spins and 2 would help greatly in clarifying the parametrisation of 
/2(1565). However, the analysis of these data must take into account contributions from a — ► Air; that has not 
been done up to the present. 

Fig. 13(c) shows the Argand loop for / (1370) in elastic scattering. The loop is cut off at the left by the 
effect of Air inelasticity. It is remarkable that the loop is very close in shape to the circle given by a Breit- 
Wigner resonance of constant width. This provides some support for the constant width approximation which is 
frequently used. The resonance is behaving to first approximation as a simple pole with appropriate widths to 
2n and An. However, Fig. 13(b) compares the phases for the s-dependent form of the amplitude (full curve) and 
constant width (dashed); these phases are measured from the origin of the Argand diagram. There is a sizable 
difference in phases, but only above the upper half-width of the resonance. 

3.11 Parametrisation of m(s) 

For convenience, an algebraic parametrisation of m(s) is given here, to allow reconstruction of the amplitudes for 
/o(1370), /o(1500) and a. It is not possible to find a simple accurate formula dealing with all three mass ranges. 
Instead, formulae will be given which are sufficiently accurate for extrapolations from the physical region to the 
poles. This implies weighting the fit to m(s) in the vicinity of these poles. 

The form of parametrisation is guided by the facts that (a) is parametrised as a Fermi function, (b) the 
dispersive term is given approximately by the gradient of this function. Then m(s) is expressed as a sum of 
terms 



Table 3 gives numerical values of parameters. Note that the a parametrisation applies only to the vicinity of the 
pole; if the full form of m(s) is needed for the a over the entire mass range, values may be read from Fig. 7 or 
the author will supply numerical tables. 

A full account of pole positions on the many possible sheets is not helpful. The sheets may be labelled by 
the sign multiplying i. Table 4 then lists a representative set for /o(1370) and /o(1500). One sees immediately 
that the sign attached to i for KK and r\r\ has little effect on the pole position, because these inelasticities are 
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Table 3: Parameters fitting m(s) in units of GeV. 



small. What matters are the signs of i for irir and 4tt sheets. The imaginary part of the pole position changes 
substantially when the sign for the 4ir sheet changes. This is a familiar effect of a strong inelastic channel. The 
experimental line width is close to the average of the results for the two 4tt sheets. 
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Tabic 4: Pole position on various sheets. 



4 Fits to elastic scattering 

Four sources of information on elastic scattering are fitted simultaneously with production data. The first are 
Cern-Munich data from irp — ► (7T7r)n. Secondly, K e 4 data of Pislak et al. are included and constrain the tttt 
S-wave phase shifts up to 382 MeV [38]. Thirdly, Caprini et al have made a prediction of nn phases using the 
Roy equations [39]. Their prediction up to 925 MeV is included with errors which are adjusted to give x 2 = I 
per point. These predictions are particularly important in constraining the scattering length and effective range. 
Amplitudes for irir isospin 2 amplitudes are also included, and parametrisations are given below. Fourthly, BES 
II data on J/ty — > cj7r + 7r~ provide an accurate parametrisation of the a pole, as discussed in Ref. [32]. The 
prediction of the lower side of the a pole from the Roy equations is precise, but the BES data determine the 
upper side more accurately because of effects arising from the sub-threshold KK and r\r\ contributions. 

In fitting Cern-Munich data, J p = + contributions are included from a and /o's at 980, 1370, 1500, 1790, 
and 2020 MeV; the last of these is above the mass range of data, which finish at 1.89 MeV, but it needs to be 
included because of its large width. The /o(1710) is dominantly ss and is not expected to contribute strongly; 
any possible contribution is absorbed into the parameters of /o(1790). 

Contributions are allowed for J p = 1" from p's at 770, 1450, 1700, 1900 and 2000 MeV, although the p(2000) 
is included only for completeness. For spin 2, / 2 (1270) plays a dominant role, but /2(1565) is definitely needed, 
as is some contribution from either or both of /2(1920) and /2(1950); the latter two however, cannot be separated 
cleanly. The 1 = 2 D-wave is also included, with formulae discussed in the next sub-section. 

For spin 3, p(1690) plays a strong role, but there is definite evidence for some additional contribution from 
,03(1990). Finally, some definite but small contribution from /4(2040) is needed. 
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Before going into details, final fits to Cern-Munich moments are shown in Figs. 14 and 15. Panels are labelled 
by L,M of spherical harmonics fitted to data. The fit is quite adequate, but the eventual x 2 is 3.13 per point. 
There is, for example, a definite systematic discrepancy with the Y(51) moment, Fig. 15. This shows structure 
around 1270 MeV which cannot reasonably be attributed to interference with the low mass tail of p(1690): the 
effect is too large. Near 1550 MeV, there are discrepancies with Y2 and Y4 moments, probably because of the 
effect of the pp threshold on / 2 (1565); this is not explicitly included. 




M (GeV) M (GeV) 

Figure 15: Fit to Cern-Munich moments with M = 1 



The fit to /o(980) is constrained within the linear sum of statistical and systematic errors quoted by BES 
for its mass, g 2 (Tnr) and g 2 (K K) / g 2 (tttt) . The last of these is a useful constraint, but the fit does optimise the 
other two parameters within the BES errors, showing there is no conflict with Cern-Munich and Crystal Barrel 
data. The mass and width of the prominent p(770) also need to be fitted freely, with the results M = 778 MeV, 
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r = 153 MeV, in satisfactory agreement with PDG values. A small detail is that its coupling to 4ir is included 
using on resonance the PDG estimate of the branching ratio; including this effect marginally improves the fit to 
the tail of the p above 1 GeV. 



4.1 1 = 2 amplitudes 

There are experimental data on ir + p — > (7r + 7r + )n and hence the 7 = 2 S-wave amplitude. Wu ct al. have 
calculated the expected inelasticity [40]. These two inputs have been fitted empirically as follows: 

f(I = 2) = ( V exp(2iS)-l)/2i, (27) 
r] = 1.0-0.5p4*(«), (28) 

tan<5 = n7+s/mi {a + bq2 + cqi + fq% (29) 

where q is pion momentum in the irir rest frame, and a = -0.0444, b=-0. 12508, c = -0.00561, d = 0.00014, all in 
units of pion masses. 

The 1=2 D-wave has been studied carefully by Pelaez [41]. His formula for the phase shift 5 is used. The 
elasticity rj is set to 1 up to 1.05 GeV, and thereafter taken from the work of Wu et al. [40] as 

»y = l - 0.2(1.0- 1.05 2 /s) 3 . (30) 
The effect of this amplitude is very small in the present analysis. 



4.2 Comments on fitted resonances 

The 7T7T S-wave plays an important role in fitting moments with L = 0, 1, 2 and 3. The value of T^-n for /o(1500) 
is well determined and significantly larger than has been generally assumed. Values of are collected into 
Table 5. Values of \ 2 are renormalised downwards by a factor 3.13 in order to allow for the high x 2 of the 
fit. The value of T2tt for /o(1790) is not well determined because of overlap with the broad /o(2020). Argand 
diagrams are shown in Fig. 16. 

Without /o(1370), a different type of solution can be obtained, but without the extra loop of Fig. 16 between 
1.15 and 1.25 GeV. If the mass and width of /o(1500) are fitted freely, the x 2 without /o(1370) is worse by 26. 
However, the width of /o(1500) goes up considerably, obviously because it is simulating missing / (1370); also 
the mass goes down. If both mass and width of / (1500) arc held fixed, x 2 is worse by 45 than with /o(1370) 
included. 

Either or both of p(1450) and /o(1700) are required to reproduce L = 1 moments, but are poorly separated. 
One can remove either with a change in x 2 < 10. The mass and width of p(1450) are taken from Babar data 
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Figure 17: The fit to Cern-Munich moments without mixing between a, /o(1370) and /o(1500) 
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Figure 18: The fit to the mr mass projection of BES 11 data for J/* — > a>7r + 7r 



[42], where a conspicuous ai(1260)7T signal is observed. The p(1900) is fixed in width to the value 145 MeV, 
the mean of the two sets of Babar data [43]. It gives a significant improvement, but is not well separated from 
p(2000). It may well be a radial excitation of p(1450), but could also be a cusp effect due to the opening of the 
pp threshold. A fit using the narrow width of Frabetti et al [44] is somewhat poorer. This raises the possibility 
of a resonance which has been attracted to the pp threshold; Frabetti et al may be sensitive to the threshold 
effect, while Babar data may be more sensitive to the resonance. 

Spin 2 states contribute strongly to moments with L = 2 to 5. The /2(1565) is definitely required but cannot 
be parametrised accurately because of missing information on . 

The p 3 (1690) makes a strong contribution to moments with L = 5 and 6. Its optimum mass is 1709 ±6 MeV. 
Reducing it to the PDG average of 1688.7 ± 2.1 docs not affect the fit to L = 6 moments significantly, but does 
make the fit to L = 5 moments worse by 20 in \ 2 ■ The value of T2n/^'tot is 0.224 ± 0.032, in close agreement 
with the PDG value. 

The fits shown in Figs. 14 and 15 include mixing between a, /o(1370) and /o(1500). Fig. 17 shows the fit 
to Y0, Y2, and YA with the mixing removed. Without this mixing, the fit to Y2 is not satisfactory, because 
interferences between / 2 (1270) and /o(1370) are not accurately fitted. Extra phase variations from the mixing 
play an important role. Such mixing is to be expected for strongly overlapping resonances. 

Fig. 18 shows the fit to the 7T7t mass projection of BES data for J/\l/ — + Luir + 7r~ , normalised to 1 at the 
highest data point. The a pole is at 470 ± 30 — i(260 ± 30) MeV. These values have changed little from the work 
of Ref. [32]. Parameters fitted to the a amplitude are shown in Table 6 in units of GeV. A parametrisation of 
m(s) in the vicinity of the a pole is shown above in Table 3. 
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State 



T 27r (McV) AF 



/ 4 (2040) 
p 3 (1990) 
Ps(1690) 
/ 2 (1565) 
/o(1500) 



16 ±7 5.3 

36 ±8 12.4 

44 ± 3 497 

46 ± 14 33 

61 ± 5 149 



Table 5: Parameters of fitted resonances; column 3 shows the (renormalised) change in x 2 when the resonance 
is removed from the fit to Cern-Munich data. 



5 Fits to pp — > r]r]7T data 

The r/Tr and r\r\ mass projections are shown in Fig. 19 for data in liquid hydrogen and gas. In (b) and (d) there 
are obvious peaks due to ao(980). In (a) and (c) there is a strong peak due to /o(1500) and a lower peak which 
is naturally attributed to /o(1370); it cannot be due to /2(1270), whose coupling to 7777 is much too weak. 

The ingredients in the fit are a (980) and a 2 (1320), a, /o(1370), /o(1500), / 2 (1270) and / 2 (1525). The 
/ 2 (1270) — > r\r\ contributes very little intensity, because of its small branching ratio, but it is desirable to include 
it in order to accomodate possible interference effects. This is done, fixing the branching ratios between 777771-° 
and 3tt° data according to the PDG value for Bi?[/ 2 (1270) 7777] /BR[f 2 (1270) tttt]. 

An earlier Crystal Barrel paper [16] has noted that there is evidence for / 2 (1525) — > 7777. However, its 
branching ratio is sufficiently low that it would not be seen in 3n° data. This result is confirmed here, and the 
/2(1525) makes a large improvement to fitting 777771 data, A% 2 = 344 summed over data in liquid and gas. 

Additional fits have been made including the exotic wi (1400) in the irr] P-wave; however, it makes no significant 
contribution. The ao(1450) is above the available mass range in 777T, but its low mass tail could contribute and 
this has been tried. However, the fits are no better than can be obtained by including form factors into the high 
mass tail of ao(980) in r/ir and KK channels. Incidentally, the decay channel ao(980) — > 77V is included, assuming 
that the ratio g 2 [ao(980) — * ?7'7r]/(j 2 [ao(980) — > 77tt] = (0.6/0.8) 2 , as predicted by the pseudoscalar mixing angle. 
This does improve the fit significantly, but needs further study for pp — ► ?77r 7r , where the a (980) signal is more 
distinct. 

The Argand diagram fitted to the 7777 S-wave is shown in Fig. 20. There is a distinct cusp between / (1370) 
and /o(1500). Masses and widths fitted to both resonances are entirely consistent with 37r° data. 

Without / (1370) in the fit, the a amplitude makes some contribution to replacing it. However, x 2 is worse 
by 317 for liquid data, i.e. > 17 standard deviations, and 68 for data in gas, > 8 standard deviations. The 
significance level in liquid is high and confirms that the peak at <~ 1300 MeV in 7777 cannot be fitted adequately 
with the contribution alone. This result by itself is sufficient to show the existence of /o(1370), but the 
determination of its mass and width are much poorer than from 37r° data. 

The fit determines the branching ratio of a — > 7777 compared with a — > tt 7t . There is, however, a substantial 
error which arises from interferences between the three 7r 7r° contributions in 37r° data. The a contribution is the 
largest one to 3tt data, and is subject to some flexibility depending on the precise s-dependence of the amplitude. 
There is a subtlety concerning the evaluation of the branching ratio. In the 3tt data, the integrated intensity 
comes not only from the three individual 7T7T combinations 12, 23 and 13, but also from interferences between 



Parameter Value 



M 0.900 

6 3.728 

61 0.092 
A 2.882 
g 4 0.012 



Table 6: Parameters fitted to the a amplitude in units of GeV 
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2 3 0.45 0.98 1.51 

Mass 2 (r]ri) (GeV 2 ) Mqss^tt) (GeV 2 ) 



Figure 19: Mass projections from data on pp — ► 777771- at rest in liquid hydrogen and gas; points show data and 
the histogram shows the fit 
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Figure 20: The Argand diagram for the r\r\ S-wave in 7/7771"° data in liquid; masses are shown in GeV. 



them. These interferences are quite large. What one needs is the branching ratio of an isolated resonance without 
these interferences. Putting this point in a different way, the resonances have specific coupling constants g 2 to 
each decay channel. One needs to derive these allowing for the fact that the three combinations interfere in the 
integrated intensity. This is done by evaluating intensities with and without these interferences. The result is 

5 2 (w)/.9 2 (7T7r) =0.19 ±0.05, (31) 

in close agreement with the earlier determination of Ref. [22], namely 0.20 ± 0.05. 

Branching ratios for /o(1500) and /o(1370) to 7777 compared with irir are subject to the same large uncertainty 
arising from the range of possible fits of a to 37r° data. Results are 

r[/ (1500) -► 7777)/r(/ (1500) -► 71-71-] = 0.135 ±0.04, (32) 
r[/o(1370)^ 7777)/r(/o(1370)^7T7r] = 0.19 ±0.07. (33) 



6 Refitting data on J/^f — > 07r + 7r 

These data provided earlier evidence for / (1370) [21]. They are refitted here using for a, /o(1370) and /o(1500) 
the formulae of Section 2 and the parameters determined here. The fit to 07T7t data is displayed in Fig. 21, 
retaining the layout of the BES publication for comparison. Data on (f>KK are refitted simultaneously, but 
changes are barely visible on figures and will therefore not be shown here. 
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M(7T + 70 (GeV/c ! ) M(tt 4 7t") (GeV/c z ) 

Figure 21: Fit to BES II data on J/^f — > 07r + 7r~. (a) Dalitz plot from data; (b) fitted Dalitz plot; (c) irir mass 
projection and fit (upper histogram); experimental background is shown by the lower histogram; (d) fit to the 
4>tt + mass projection with the same format as (c); (e) the coherent sum of + contributions (full histogram); the 
lower histogram shows the a component; (f) the contribution from /2(1270). 



Technical details are as follows. It is important to constrain the branching ratios of /2(1270) and /o(15QO) 
between KK and tttt within one standard deviation of PDG averages, so as to stabilise the fit to <pKK data. 
The new a parametrisation replaces the previous one, but has almost no effect: a 0.7 standard deviation change 
in log likelihood; there are much larger changes from reparametrisation of /o(1500) and /o(1500). The a — > KK 
contribution is very small and has little effect on the fit. 

A reminder is necessary concerning the 4>tttt data. The BES analysis located a peak in <fiir at 1500 MeV, 
arising from the triangle graph due to J/ VP — > K + K~tt + tt~ , followed by rescattering KK — > <fi. This was 
eliminated by a kinematic cut, creating the dips in the 07r mass spectrum of Fig. 21(d) at 1.5 and 2.8 GeV. 

There is a peak at 1350 MeV in the tttt mass projection (c) attributed to interference between /2(1270), 
/o(1370) and /o(1500). The new fit is shown by the upper histogram; the lower one shows experimental back- 
ground. The fit is marginally poorer than that of the BES publication, which had the freedom to optimise the 
mass of / (1370) to 1350 ± 50 MeV and the width to 265 ± 40 MeV. It resulted in an /a (1370) contribution 
larger than /o(1500), though both were needed. Their roles are now reversed, with a 4.6% intensity contribution 
to </>7r + 7r~ from /o(1370) and a 6.1% contribution from /o(1500). On Fig. 21(e), the full histogram shows the 
coherent sum of J p = + contributions and the lower dashed histogram the a contribution. Fig. 21(f) shows 
the /2(1270) contribution. It is cleanly separated from + by angular distributions for production and decay of 
resonances. 
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M(7T*7T) (GeV/c 2 ) M(nV) (GeV/c 2 ) 

Figure 22: (a) Fit to <^7T7r data without /o(1370); (b) the fit with a very broad + signal at 1650 MeV interfering 
destructively with / (1500). 



Fig. 22(a) shows the fit without /o(1370). The narrow /o(1500) is unable to compensate via interferences 
with er, /2(1270) and /o(980) and the resulting fit is visibly poor. The fit is worse than Fig. 21 by 34.8 in log 
likelihood, i.e. > 8 standard deviations. This is sufficient to confirm the presence of /o(1370), but not enough 
to influence its fitted mass and width. 

There is evidence that the /o(1370) is resonant. If its phase variation is artificially suppressed, leaving the 
line-shape unchanged (even though the resulting amplitude is non-analytic, therefore illogical), the fit is worse 
by 12.4 in log likelihood, ~ 5 standard deviations. 

It has been suggested that the /o(1790) in BES data might be eliminated by fitting a broad J p = + nir 
signal with which /o(1500) interferes destructively, producing an interference minimum at ~ 1600 MeV. This 
suggestion was based on the possibility that the high mass tail of the a might contribute strongly. That now 
appears illogical, since the a amplitude falls rapidly with increasing s and is already small at 1 GeV. 

The suggestion has been tested by fitting with a broad component with a width 1000 MeV and an optimised 
mass of 1650 MeV. Fig. 22(b) shows the poor fit. It is of course possible to tune the width of this broad 
component to a lower value and produce a reasonable fit. However, it makes little sense to invent a new broad 
component when there is independent evidence for / (1790) with consistent parameters in J/vt — * j(Att) data 
from both Mark III [45] and BES I [46]. Data on J/* -> 7^ [47] exhibit a striking 0+ peak at 1812 MeV, 
which is consistent within errors with the upper side of /o(1790); the u>4> threshold is at 1801 MeV. 

7 Other data 

Let us recall that the first evidence for /o(1370) came from data at the Argonne and BNL laboratories on 
7T7T — > KK in the 1970 era. There is a distinct dip between the 1300 MeV peak and the narrow /o(980). 
Fig. 10(a) of Ref. [22] shows that without /o(1370) this dip cannot be fitted: the /o(980) is too narrow and 
interferences with the broad a and /o(1500) fail to fit the data. Though these data alone at not sufficient proof 
of the existence of /o(1370), they are entirely consistent with parameters fitted here. The fit made in Ref. [22] 
has been rerun using the new parametrisations of the a, /o(1370) and /o(1500) determined here. There is no 
significant change to conclusions of Ref. [22] . The analysis of Crystal Barrel data on KKtt channels by Anisovich 
and Sarantsev [30] finds /o(1370) mass, width and coupling constant in acceptable agreement with those reported 
here. 

Data of Barberis et al. [48] on central production of 7T7t reveal quite different azimuthal distributions for 
/o(1370) and / (1500), suggesting that two separate resonances are needed, whatever the nature of / (1370), 

8 The need for further analyses 

Those who question the existence of /o(1370) should be concerned also about the existence of ao(1450). Despite 
its appearance in the Summary Table of the Particle Data Book, it is subject to the same questions about 
dispersive corrections as /o(1370). A fresh analysis is needed of data on pp — > r)Tr°ir°, where it was discovered 
[49], including m(s) into the parametrisation. It has also been observed in pp at rest — > u/pn in the top channel 
[50]. The threshold for this final state will contribute strongly to m(s), as may the possible decay channel 
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a (980)cr. A combined analysis of data on 777777° and ojpw is in progress and will be reported separately. If the 
ao(1450) survives, it is very likely to be a qq state, since none other is available for the isospin 1 component of 
the nonet in this mass range. If so, it is plausible that /o(1370) is likewise qq (mixed with the glueball). 

Next, it is highly desirable to fit all existing data for the J PC — 1 sector including dispersive corrections. 
A start has been made on this type of analysis by Weng et al. [51]. The greatest need is for data on tttt — > 4tt. 
These are needed to pin down details of the 4tt final state with J PC = ++ , 1 and 2 ++ . Such data were 
presented in a preliminary form by Ryabchikov at Hadron95, but have not been the subject of a full publication 
yet. The analysis including full dispersive effects is doubtless a major undertaking, but even limited information 
about the 47r channel would be very important. 

Yet another example which may be affected strongly by dispersive effects is the J — 1~ + sector. There 
is substantial evidence for tt 1 (1600) in decays to &i(1235)tt, t/tt, / 1 (1285)rr and pir. There is also evidence for 
structure in r/ir at <~ 1400 MeV. This is close to the sharp thresholds for & 1 (1235)7r and /i(1285)7r. It is still an 
open question whether there really is a resonant 7Ti(1400), or whether it can be fitted adequately as a threshold 
effect. It is also possible that there is a resonance associated directly with these thresholds. What is needed is 
an analysis including dispersive effects due to the thresholds. 

9 Concluding Remarks 

The /o(1370) is highly significant statistically in 5 sets of data: pp — > 37r° at rest in liquid hydrogen and gas, 
corresponding data for the 777771"° data channel, and J/^> — > </>7r + 7r~. Overall, it is statistically more than a 52 
standard deviation effect. What is also important is that fitted parameters of /o(1370) agree remarkably closely 
between the 37r° data in liquid and gas. There is weaker but consistent evidence for / (1370) in 7T7t — ► KK. 

The data cannot be explained by the high mass tail of the a, because it is too broad. The requirement 
for a peak in data with a full width of 207 MeV requires an additional narrower state, identified here with the 
/o(1370). No such pole has appeared within the present parametrisation of a — > 47r. 

Dispersive contributions due to the opening of the 47r channel have large effects, renormalising the Breit- 
Wigner denominator. This severely limits the range of r4 7r /r 27r which can be successfully fitted to data. Despite 
the strong effect of the 47r threshold, the resonance loop on the Argand diagram is very close to a circle. This 
provides some justification for the common practice of fitting a Breit-Wigner resonance of constant width: the 
resonance behaves to first approximation like a simple pole. However, the phase is significantly affected once 
one reaches a mass more than one half-width above resonance. A simple Breit-Wigner amplitude is adequate for 
finding resonances; including m(s) uncovers the dynamical effect of the threshold. 

There are presently no significant inconsistencies in parameters of /o(1370) between sets of data on 7T7r, KK 
and 7777 channels. Wide variations of mass and width appear only in analyses of 47r data. However, those analyses 
do not presently include a — > 4rr, which is found here to play a large role. Unfortunately, these analyses need to 
be repeated including the a contribution. Disagreements therefore exist only about schemes into which different 
authors wish to fit the known states. Those schemes should not be used as the basis for claiming that /o(1370) 
or any other resonance does not exist. 

Some authors, for example Maiani et al. [29], raise concerns about mass differences between / (1370), 
7X0(1430) and ao(1450). It now appears that /o(1370) is nearly degenerate in mass with ai(1260), /i(1285) 
and /2(1270). So it appears to pose no particular problem. It is likely that the mass of yYo(1430) is influenced 
by its strong coupling to Krf, whose threshold opens at ~ 1450 MeV. It is also known that a (1450) appears 
only weakly in the rjir channel; it is likely that its mass is affected strongly by coupling to up and ao(1450)cr 
thresholds, just as the mass of /2(1565) is close to the uu threshold, and much lower than the mass of a2(1700). 

The present analysis of elastic data produces new quantitative estimates of r 27r for /o(1500), /2(1565), 
P3(1690), p3(1990) and /*4(2040). Data on 7T7t — ► 47r would help greatly in confirming the present analysis, and 
parametrising more accurately 1 states and / 2 (1565). 

A final speculative remark emerges from the mixing between a, /o(1370) and /o(1500), which appears to 
be a necessity in fitting the elastic 7T7t data. This mixing suggests a possible analogy with chemical binding. 
In the hydrogen molecule, two configurations of protons and electrons mix. This is the familiar process of 
hybridisation. The lowest eigenstate may be calculated using a Variational Principle. The suggestion made here 
is that confinement involves a hybridisation due to overlapping of nearby resonant states. One linear combination 
of states with favourable SU(3) configuration moves down and other combinations are pushed up in energy, in 
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a way analogous to formation of the covalent chemical bond. A similar Variational principle is involved in the 
formation of a superconductor. 

The relevance of such an idea to the confinement process itself is a matter of conjecture without calculations 
to support it. The idea is that there is feedback between the formation of resonances and a dynamic confinement, 
i.e. condensation. Such a mechanism would explain naturally why most of the huge number of possible molecular 
states are not observed: most of them are driven upwards and become a continuum. 
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